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Abstract. One of the central open problems in nuclear physics is the construction 
of effective interactions suitable for many-body calculations. We discuss a recently 
developed approach to this problem, where one starts with an effective field theory 
containing only fermion fields and formulated directly in a no-core shell-model space. 
We present applications to light nuclei and to systems of a few atoms in a harmonic- 
oscillator trap. Future applications and extensions, as well as challenges, are also 
considered. 
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1. Introduction 

How does the complexity of nuclear physics arise from the relatively simple QCD 
Lagrangian? Of the open questions in nuclear structure, this is one of the broadest 
and most fundamental. The long-range goal of nuclear-structure theory is indeed to 
calculate the properties of finite nuclei starting from the strong-interaction physics of 
QCD. Although this goal is still several years away, a confluence of nearly parallel 
conceptual and computational developments is stirring unprecedented excitement about 
our ability to reach it. 

Tremendous progress has been made in the last decade P, [2] in our understanding 
of how nuclear structure arises from the properties of the interactions among nucleons 
inside a nucleus, even though these interactions are often modeled only in terms of 
ad hoc potentials. In the previous twenty years a number of many-body techniques 
has been developed for exactly solving the nuclear few-body problem [3l H]. The 
application of these approaches to nuclei ranging in mass from A = 2 to A = 15 
has clearly demonstrated how the structure of these light nuclei, i.e., binding energies, 
excitation spectra, electromagnetic moments, etc., arise directly from the properties 
of nucleon-nucleon (NN), smaller three- nucleon (NNN), and perhaps some tiny four- 
nucleon (NNNN) interactions. One of the most tantalizing current problems is how to 
extend the many-body approaches for light nuclei to medium- and heavy-mass systems. 

Independently, a framework has been developed based on effective field theories 
P, E] to construct nuclear potentials that respect the symmetry pattern of QCD and 
produce observables in a systematic and controlled expansion in powers of momentum. 
These ingredients restrict the shape of nuclear interactions, particularly in the range of 
pion exchange, and encapsulate the complicated short-range physics into a number of 
"low-energy constants". Now, the first results are emerging [7], where the low-energy 
constants, so far simply fitted to data, can be calculated in full lattice QCD simulations. 
How can one link these results to nuclear properties? 

In this paper we map some of the landmarks and crucial crossings in the long road 
from QCD to nuclear structure. 

2. Effective Interactions and Operators 

The general goal of microscopic nuclear-structure theory is to begin with the free inter- 
nucleon interaction and determine, using many-body quantum mechanics, the properties 
of finite nuclei. In principle, one would like to solve the many-body Schrodinger equation 
for all A nucleons, 

=K|*a), (1) 

in the "full" Hilbert space S. Here the Hamiltonian is 

A A A 

^ = E^i + E% + E V^Jk + ..., (2) 
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where the first term involves the kinetic energies ti and the subsequent terms, the 
potentials among an increasing number of nucleons. It used to be assumed that only 
accurate NN interactions were required, although evidence has now accumulated for the 
need for NNN interactions [8] (and perhaps NNNN interactions). Relativistic effects do 
not play a significant role in understanding low-energy nuclear structure, and can be 
included as corrections around the non-relativistic limit. 

In general, Eq. ([1]) cannot be solved in the full Hilbert space S, because of the 
infinite number of configurations, so the problem must be truncated to a smaller Hilbert 
space 5" of dimension d (called the shell-model basis space or simply the model space). 
Because the full space has been truncated, full space operators cannot be used, but must 
be replaced with effective operators appropriate for the given size of the model space. 
In this case, H must be replace by the effective Hamiltonian H', such that 



where = P|\E'^), P is a projection operator from S into S', and the Ef^ are a subset 
of dimension d of the exact eigenvalues in Eq. ([T]). The projections $^ are usually not 
orthogonal, so one must construct the biorthogonals such that 



Using the biorthogonals, one can easily obtain the H' that satisfies Eq. ([3]), i.e., 



It should be noted that H' will usually be non-Hermitian, because of the non- 
orthogonality of the l^js). All other physical operators relevant to the nuclear system 
being investigated, e.g., the rms radius operator, the electromagnetic-moment operators, 
the transition operators, etc., must be renormalized in a similar manner for use in the 
given model space. 

In heavy nuclei, a drastic projection is needed in the standard nuclear shell model 
(SNSM), consisting of an inert, closed-shell core and a few valence nucleons. The 
theoretical construction of SNSM effective interactions and operators, both empirically 
and microscopically, has an extremely long and large history in nuclear physics [9], [10] . 
The microscopic approach has been more-or-less unsuccessful due to problems connected 
with the convergence of the perturbation-theory expansion for the effective operators 
and because of the complexity of the calculations. 

A response to these difficulties was found in the no-core shell model (NCSM) 
[H [21 [H] , in which all A nucleons in the nucleus are active. In its usual formulation, 
the NCSM employs a harmonic-oscillator (HO) single-particle basis of frequency u, 
and nucleons are allowed to share a maximum number N^ax of oscillator quanta (the 
maximum principal quantum number of the wave functions) above the minimum-energy 
configurations. It is important to include all configurations up to some total energy. 
This guarantees that all spurious center-of-mass components will be projected from 
the final wave functions, when using the Lawson procedure [12]. The need for effective 



H'\^0) = Epl^fs), 



(3) 




(4) 



H' = Y.\<^p)Ep{^p 



(5) 



Effective interactions for light nuclei 



4 



interactions and operators, in order to calculate nuclear properties, exists whether one is 
performing SNSM or NCSM calculations. However, the physical nature of the required 
effective operators will be quite different in the two cases. In the simpler NCSM, it 
involves only unitary transformations based on Eq. ([5]) [TTl [T3] . 

Although the determination of the effective operators for NCSM calculations is more 
straightforward, the generated effective operators are still dependent upon the nature of 
the inter- nucleon interactions employed in the calculations [H] , a low-body (usually only 
two-body) cluster approximation, and the size of the model space. Thus, it is highly 
desirable to find a better method for calculating these NCSM effective operators, which 
is not potential dependent, yet is directly related to the underlying QCD symmetries 
and can always be employed in small model spaces. 

3. Effective Field Theory and the NCSM 

The framework of effective field theory (EFT) emerged in the late 1970s [15] from the 
realization that, because in a quantum field theory any physics problem is a many-body 
problem, one never has access to the "full" Hilbert space 5*. It is always necessary to 
truncate the Hilbert space so as to exclude states associated with energies beyond those 
we can confidently access experimentally and theoretically. This paradigm shift implies 
that interactions are only defined in the context of a model space: the Platonic concept 
of "the" Hamiltonian is not particularly useful in physics. 

What is needed is a method to construct effective interactions and operators 
whether or not the underlying, higher-energy physics is known. Since relativistic 
quantum effects of arbitrary complexity exist unless they are forbidden by a symmetry, 
any effective interaction or operator with assumed symmetries should be included. Still, 
the intrinsic connection between Hamiltonian and model space would be problematic 
without a way to ensure that the arbitrariness in the choice of model space does not 
contaminate observables. The projector P always contains at least one dimensionful 
parameter, the ultraviolet (UV) cutoff A, for definiteness taken here to be a momentum. 
One thus requires that the Hamiltonian H' depend on A in such a way that observables 
at momenta Q <C A are independent of how P is chosen, and in particular, independent 
of A. This is termed renormalization-group (RG) invariance. 

Nuclear physics is a great arena for these ideas, because it is not easy to find 
solutions for QCD at momenta Q < Mqcd ~ 1 GeV. Yet we know the symmetries 
of the QCD Lagrangian. In particular, chiral symmetry plays an important role, even 
at low energies, thanks to the appearance of light pions, due to spontaneous and small 
explicit breaking. It is, thus, not difficult to write the most general Hamiltonian with the 
appropriate degrees of freedom and symmetries [5] , which is a generalization to systems 
with more than one nucleon of the Lagrangian used in chiral perturbation theory. 
Interactions among nucleons consist of pion exchanges and contact interactions, which 
subsume short-range dynamics (say, exchange of heavier mesons). At very low energies 
even pion exchange can be treated as short ranged, leaving only contact interactions in 
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the theory. 

The real challenge is ordering the infinite number of interactions so that observables 
can be calculated in an expansion in powers of Q/Mqcd- A guide is provided by the RG: 
a truncation of this expansion at any given order must respect RG invariance except 
for small errors contained in higher orders. This approach had mostly been applied 
in particle physics to systems where unitarity could be accounted for perturbatively. 
In nuclear physics, the leading order (LO) must contain non-perturbative physics to 
generate nuclear bound states and resonances. Only subleading-order corrections, if 
truly corrections, should be treated in perturbation theory. 

Most applications of EFT in nuclear physics have been carried out using a 
continuum free-particle basis [6]. This facilitates the study of reactions, in addition 
to structure, but unfortunately, with foreseeable computational resources, is limited to 
A < 4. To increase A, we need to limit the number of accessible one-body states 
by introducing an additional, infrared (IR) momentum cutoff A, which discretizes 
momentum. This is a traditional method in QCD itself, where simulations are carried 
out on lattices, and in addition to an UV cutoff A ~ 1/a, with a the lattice spacing, 
there is also an IR cutoff A ~ with L the lattice size. Lattice regularization can 
be, and has been, applied to nuclear EFT as well [161 [13 HH] > where it is particularly 
suited to the study of nuclear matter at finite temperature. 

The successes of the SNSM suggest, however, that formulating the nuclear EFT 
in an HO basis might be an efficient way to reach larger, finite nuclei. An UV cutoff 
equivalent to the cutoff used in a free-particle basis can be defined in an HO basis: if fi is 
the two-body reduced mass, A = '^2fi{Nmax + S/2)uj is the momentum associated with 
the energy of the highest HO shell included explicitly in the two-body system in the 
center-of-mass frame. The HO frequency cu defines an additional IR cutoff A = V^/b, 
where 6=1/ ^/JIuJ is the HO length, which plays a role analogous to the box size L in a 
lattice discretization [17]. 

EFT in an HO basis is essentially the NCSM formulated directly within the model 
space: EFT provides the form of the effective interactions and operators needed in 
the NCSM. Since NCSM is a full diagonalization approach in a basis constructed with 
HO wave functions, each model space is determined by two parameters: Nmax and u, or 
alternatively, A and A. Consequently, all low-energy constants are functions of these two 
parameters. As in any EFT, they need to be determined either from some experimental 
data or, eventually, from QCD itself. Enormous progress in lattice QCD has already 
produced [7] NN scattering lengths, albeit at unphysical values of the pion mass. It is 
reasonable to expect that the next few years will see scattering lengths at lower pion 
masses, effective ranges, other NN quantities, and perhaps even observables involving 
more nucleons. A key aspect of this work is the use of Liischer's formula [19], where the 
energy levels inside a box are linked to scattering parameters. 

The drawback of formulating the EFT in an HO basis is that the connection to 
scattering states becomes less obvious than when using a continuum free-particle basis. 
As we show in the rest of this paper, this is not an impassable roadblock, for two reasons. 
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First, one can always determine the low-energy constants from bound-state data 
{e.g., binding energies). We illustrate the method in Sec. HI where we use data from 
a few bound states to fix the low-energy constants and then many-body theory in the 
form of the NCSM to calculate properties of larger systems [20]. It is true, however, 
that as the EFT order increases, the low-energy constants multiply, and it is desirable 
to determine them using the more abundant scattering (experimental or lattice) data. 

Second, it is possible to connect the energies inside an HO well to scattering 
parameters in a way [21] very similar to Liischer's formula. We show how this can 
be done in Sec. |5] in the simpler problem of spin-1/2 fermions without isospin [22], 
which is relevant for trapped two-state atoms close to a Feshbach resonance. 

4. Light Nuclei 

The first direct application of EFT principles to the derivation of effective interactions 
in finite NCSM model spaces was presented in Ref. [20]. In that work, we have opted 
for a theory without explicit pions for two reasons: (i) for very-low-energy processes 
(involving momenta less that the pion mass) the formalism becomes very simple, and 
(ii) the same techniques would be readily applicable to the pionful EFT. 

In pionless EFT at LO, the Hamiltonian H' in Eq. ([3]) can be written as a sum 
of the relative kinetic energy, two contact interactions in the ^5*1 and ^5*0 NN channels, 
with corresponding parameters Cg and Cq, and a contact three-body interaction in the 
^5*1/2 NNN channel, with parameter Dq: 

H' = j^T. (pi - p^f + Co E m - r,) + c° E '^(-1 - 

[i<3\ [i<jV [i<jf 

+ Do m-rj)6{fj-n), (6) 

[i<j<k] 

where [i < j] denotes all pairs of particles, [i < JY pairs in the S'-wave NN channel 
of spin s and [i < j < k] triplets in the spin-1/2 S'-wave NNN channel. Unlike in 
the conventional NCSM approach based on a unitary transformation, in this approach 
the interaction in each model space has the same structure, i.e., matrix elements of 
the contact two- and three-body interactions. Only the coupling constants differ in 
each model space; they depend on Nmax and lu: CQ^Nmax,^^), Do^Nmax,^)- The three 
coupling constants are fixed for each model space so that three observables, the deuteron, 
triton, and ^He binding energies, are simultaneously reproduced. With the Hamiltonian, 
thus, defined, we have investigated the energy of the first (0+, 0) excited state of ^He as 
well as the ^Li ground-state energy. 

Because the errors associated with the terms neglected in the LO Hamiltonian 
H' decrease with increasing the ultraviolet cutoff, large values of A are desirable. 
While increasing N^ax niakes the calculation impractical even for a small number of 
particles, increasing uj is not a good option either, as this would increase the errors 
associated with the infrared cutoff. Instead, for fixed a;, we have calculated the energies 
of interest in the four- and six-nucleon systems for the maximum Nmax we were able 
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to handle, extrapolating to other values of A using Eq{uj) + A{lj)/A, an expression 
motivated by a similar running in the two-body system in the continuum [6]. Then for 
fixed A, we extrapolate energies to — 0, which eliminates infrared errors, using a 
quadratic formula. Using this approach, we have predicted that the energy of the first 
(0"'',0) excited state of ^He is 18.5 MeV, in very good agreement (within 10%) with 
the experimental value. (The agreement can be understood if we consider that this 
state is very close to the four-nucleon continuum threshold, a regime well within the 
limits of applicability of the theory.) A similar analysis for the ^Li ground state also 
produced results within 30% of the experimental level, with an underbinding of about 
9 MeV. Better results can be obtained if one includes an extra term to the running, 
i. e., Eq{uj) + A{uj) /A + log(-B(a;)A) / A. While the extra term essentially does not change 
the result in the four-body system, for ®Li the new results overbind the system by only 
15%. In Fig. [1], we present the results for ^Li; they include the log term in the fit. 
Obviously, more investigations are necessary to pin down the running of the observables 
in the many-body system. 

5. Few- Atom Systems in Traps 

Beyond LO (or even at LO in the theory with pions [23]), the procedure discussed in the 
previous section to determine low-energy constants becomes impractical: the number 
of required input observables increases rapidly, and an ever larger number of states in 
light nuclei is needed to adjust the interaction. This motivated us to devise another 
approach, in which the two-body renormalization is realized at the two-body level |22j . 
In such an approach, the low-energy scattering properties in free space are related to 
the energy spectrum of two interacting particles in a harmonic trap. 

The perfect testing grounds for such an approach are systems of two-component 
fermions with a large two-body scattering length 02 in an external harmonic trap, which 
can be realized experimentally with atoms trapped by lasers in variable magnetic fields 
[23]. In this case, there is no need to take the limit a; — > 0, since u, or equivalently 
the trap length b, is given by the trapping laser. In the unitary limit 6/02 — 0, the 
spectrum for three particles is known [25] . 

In this case, the Hamiltonian is given in LO by 

+ 2/i62Co E ^i^i - rj) \ . (7) 

[i<j] J 

There is just one two-body coupling constant Co{Nmax, ^) and three-body forces appear 
only at high orders. 

In order to renormalize the two-body interaction in a finite model space, we first 
consider the two-body system in the center-of-mass frame, where the relative position 
is denoted r. Since the contact interaction connects S states only, higher angular- 
momentum states are undisturbed from HO ones. It is sufficient to consider the two- 
body wave function described by a superposition of all 5* states with the HO quantum 
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Figure 1. Ground-state energy Egs of ^Li as a function of tlic UV cutoff A (top 
panel) and of the HO frequency uj (lower panel), all in MeV. In the top panel we 
show calculated values at various frequencies (indicated in the legend) and fits of the 
form discussed in the text. In the lower panel we show values at various UV cutoffs 
(indicated in the legend) and fits of a quadratic form. 



number N = 2n < Nmaxi 

N max 1 2 

^(f)= ^n0n(r-), (8) 

n=0 

with the S'-wave HO state with radial quantum number n and An a set of complex 
coefficients. The Nmax/'^ + 1 unknown coefficients An and the energy spectrum can be 
found by solving the Schrodinger equation for the relative motion, once Cq (-^max ; '-'-') IS 
determined. 

The coupling CQ^Nmaxi^) can be determined by the following procedure. In order 
for the Schrodinger equation with the delta function to be well defined, Co(A^max, ^) and 
the energies e{Nmaxi^) (in units of u) have to satisfy an RG condition [22], 
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Figure 2. First excited-state energy of tire two-body system, e, in units of cj, as 
a function of the UV cutoff A, in units of the inverse of the trap length b. Three 
interaction strengths determined by the indicated values of 6/02 are shown. The 
ground-state energy given by Eq. ()10p is used in each case to fix the coupling constant. 
The exact values for the first excited-state energy are displayed with dotted lines. 



where L^^^\0) is a generalized Laguerre polynomial at the origin. In the A — 00 limit, 
one can show [26] that this condition leads to 

r(3/4-e(oo,a;)/2) ^ _^ 

r(l/4-£(cx),u;)/2) 202' 
That is, the energy spectrum inside the trap depends only upon the scattering length 
in units of the trap length. This formula, the HO oscillator counterpart of the Liischer 
formula for a box with sharp boundaries, was first derived in Ref. [21], using a 
pseudopotential. The CQ^Nmax,^^) can then be fixed so that in each model space one of 
the states obtained by the exact diagonalization of the Schrodinger equation reproduces 
the corresponding value given by Eq. ( fTOj) . For simplicity, we match the lowest state. 
Because just one coupling constant needs to be fixed in LO, all the other levels can be 
calculated. They satisfy Eq. ([9]) and deviate from the exact value given by Eq. (fTOll . but 
the error decreases with increasing the size of the model space. As an illustration, in Fig. 
[2]we present the running of the first excited state of the two-body system as a function of 
the dimensionless quantity Ab for three selected values of the b/a2 ratio. In all cases, the 
energy goes to the exact value with errors that decrease as 1/Ab. Faster convergence 
can be achieved by considering terms beyond leading order, that is, corrections that 
involve derivatives of the contact interaction and account for the effective range, the 
shape parameter, etc., introduced by the model-space truncation [26] . 

With the two-body interaction fixed in a truncated two-body model space, we 
turn now to the few-body problem, considering the unitary regime as well as a general 
non- vanishing 6/02 value, for both positive and negative scattering lengths. The three- 
body solutions are obtained by a diagonalization in a finite model space, constructed 
as anti-symmetrized three-body states of HO wave functions (for details on the basis 
construction in each model space, see Refs. [201 ES]). In Fig. [3l we present the running 
of the energy of the lowest three-body state (which has orbital angular momentum 
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Figure 3. Ground-state energy of the three-body system in units of w as a function of 
the UV cutoff A, in units of the inverse of the trap length b. Three interaction strengths 
determined by the indicated values of 6/02 are shown. The two-body ground-state 
energy given by Eq. (jlOp is used in each case to fix the coupling constant. The exact 
value at unitarity [25] is displayed with a dotted line. 



L = 1), at unitarity and at 6/02 = ±1. At unitarity, where a semi-analytical result 
exists [25], our results are within a few percent of the exact results even for relatively 
modest model-space sizes. However, it is worth noting that the running is rather slow, 
having in general a l/{Ab)" dependence, with a a power that is state dependent. For the 
ground state at unitarity, we have shown by direct fit that a = 1 [22], and in the limit 
Ab — i> 00, we obtain the exact result [25]. Faster convergence, although not necessarily 
faster running (that is, same power a, but larger coefficients in front of 1/(A6)°'), can 
be achieved if one introduces higher-order terms [27] . 

In the non-interacting limit {b/a2 —00), the lowest state of three two-component 
fermions has L = 1 (negative parity). As b/a2 increases, this state persists as the 
ground state beyond unitarity. However, around b/a2 ~ 1.5, the lowest state becomes 
an L = state [22], which remains lowest as b/a2 increases further. In the limit of large 
b/a2, which approaches the untrapped case, the energy of the three-body ground state 
becomes ^ — l/2yua2- This result suggests that the system of three fermions is near 
the threshold of the S'-wave scattering of one fermion on the bound state of the other 
two. 

Further investigation has shown that the known results in the four-body system [28] 
are also reproduced, although with less accuracy because of the complexity of solving 
the four-body problem in larger model spaces. Nevertheless, the errors can be reduced 
by introducing higher-order corrections, and preliminary calculations show very nicely 
this trend [221 [27]. 
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6. Future Applications and Challenges 

As discussed in Sec. 3, our goal is to obtain effective interactions and operators in 
truncated model spaces for solving the few-body problem (for fermions and bosons), 
based on the underlying assumed QCD symmetries, thus bypassing, for example, 
the need for a phenomenological potential. As we have shown, this approach works 
reasonably well for light nuclei (Sec. 4) and even better for few-fermion systems in an 
HO trap (Sec. 5). 

Now that we have achieved good results linking few-fermion energies in an HO 
trap to scattering parameters, we can pursue our original idea of applying this method 
to nuclei. Doing this is, indeed, a great challenge, as we have indicated in Sees. 3 
and 4: first, because the number of low-energy constants increases significantly as the 
EFT order increases; and, second, because the many-body calculations {i.e., within the 
NCSM) become more and more difficult as the model-space size, as defined by N^ax, 
increases for larger values of A. For given computational resources, it is, thus, important 
to obtain better converged results in smaller model spaces for NCSM calculations of 
heavier nuclei. 

Obviously, we want to extend this approach to higher orders and include the 
effective range. In preliminary calculations for few-fermion systems in a trap [261 1221 we 
have included corrections to the potential up to N^LO and observed an acceleration of 
the convergence for energies. Corrections beyond LO are treated as perturbations and 
these preliminary results show an excellent agreement at unitarity with known results 
for the three-fermion system [25] and with other methods for the four-fermion system 
[28] . One finds that this perturbative treatment of subleading interactions, demanded 
by RG consistency, gives faster convergence than a non-perturbative treatment. 

Another new idea deals with making sure that the available energy in the many- 
body system is larger than the maximum two-body energy employed in determining 
the two-body interactions. This is achieved by introducing a different total number 
of oscillator quanta A^^i each, n-body space. Preliminary calculations [27] for the 
problem of few-fermion systems in an HO trap indicate that one obtains significant gains 
in the converged results in smaller model spaces, when for a fixed value of N^l^, the 
many-body model space is increased until convergence. The next step is to apply this 
approach to few nucleon systems. 

Although we have so far concentrated on binding energies, other bound-state 
observables can be calculated with our method. Using similar techniques, we can 
construct other operators for describing observables of interest. It will be instructive to 
compare with results obtained in the traditional NCSM |14j . 

Finally, an important long-term goal of this program is to eventually complete the 
link from QCD to nuclear-structure observables. In the direction of QCD, we should 
determine the low-energy constants of the EFT expansion from lattice simulations. In 
the direction of heavier nuclei, we need to understand the limits of the pionless theory, 
and presumably go beyond it. The pionless EFT is valid for low momenta (Q < m^). 
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and so is Eq. (fTUl) . One of the challenges of future applications to nuclear systems is to 
extend the equation that determines the spectra of two trapped particles for momenta 
large enough to discern pions. 

7. Summary and Conclusions 

We have presented a discussion of the first implementation of EFT principles directly 
into a many-body method, the NCSM. We have proposed that using EFT to construct 
effective interactions in restricted model spaces used in NCSM calculations might provide 
an important step in the long road from QCD to nuclear structure. 

In one implementation, we have determined the low-energy constants by a direct fit 
to ground-state binding energies in two-, three- and four-nucleon systems, and predicted 
an excited state in ^He, as well as the ground-state energy of ^Li [20]. We have presented 
an improved extrapolation for the latter in Sec. HI For both states, results are within the 
expected errors of the pionless EFT. Such an application becomes quickly impractical 
if one considers the subleading orders in the pionless EFT, or the pionful theory, as the 
number of low-energy constants that have to be determined by bound states increases 
considerably. 

In a second implementation, we have considered a system of trapped fermions |22j . 
Such an approach would allow the determination of low-energy constants in the two- 
body system alone, as shown in Sec. [5l We have illustrated the application of the 
two-body renormalization by computing the spectrum of three spin-1/2 fermions in a 
trap at unitarity, where a semi-analytical calculation exists [25]. Thus, we have shown 
that in LO our results converge to the exact ones, in the limit of large UV cutoffs. 
The advantage of our approach is that it can be extended to finite two-body scattering 
lengths, the convergence of which we examined here, and to include other scattering 
parameters [26| [27]. 

Despite the different underlying physics, the systems of trapped cold atoms near 
Feshbach resonances and of nucleons at low energies are quite similar. One can hope 
that the same procedures can be transferred to the nuclear many-body problem, as 
we discussed in Sec. |6l thus providing a QCD-based solution to the open problem of 
constructing nuclear effective interactions. Work in this direction is in progress. 
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